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The incommensurate phase of a bilayer quantum Hall state 
is found to have a "rippled" dipole charge density whenever 
the layers are unbalanced. This tunable dipole-density-wave 
instability could be detected by sensitive capacitance mea- 
surements and by anisotropic transport. We demonstrate this 
explicitly by carrying out a Hartree-Fock calculation of the 
layer densities and capacitance for a double-layer quantum 
Hall state at a total filling factor of 1. 

PACS numbers: 73.43. Cd, 64.70.Rh, 71.10.Pm, 71.45. Gm 



I. INTRODUCTION 

The combination of reduced dimensionality and strong 
interparticle interactions can have spectacular effects on 
the nature of the ground-state and dynamical proper- 
ties of many-particle systems. This is espteeially evi- 
dent in the fractional quantum Hall regime,™! where a 
strong magnetic field is applied perpendicular to a two- 
dimensional electron gas at very low temperatures. The 
powerful magnetic field quenches the kinetic energy of the 
electrons, so that interactions between electrons domi- 
nate the energetics. The result is a highly correlated, 
incompressible quantum-liquid ground, state that sup- 
ports fractionally charged excitations.™ Even at integer 
filling factors (especially i*r=l), the combination of the 
quenched kinetic energy of interacting electrons plus ex- 
tra electronic degrees of freedom (i.e., spin or multiple 
layers) can give rise to correlated ground states that sup- 
port remaEkable topological excitations, such as charged 
skyrmions.QH 

Another important class of systems in which re- 
duced dimensionality and interactions strongly affect the 
ground state and dynamics are systems with charge- 
density-wave or spin-density-wave ground statesp and 
also systems that, exhibit commensurate-incommensurate 
(CI) transitions.LrH Systems with charge- or spin-density- 
wave ground states occur most famously (although not 
exclusively) in quasi-one-dimensional materials, where 
the reduced dimensionality enhances fluctuation effects 
and leads to broken-symmetry ground statesE Systems 
exhibiting CI transitions have broken translational sym- 
metry states with rich structures, most notably arrays 
of domain walls that arise from the competition be- 
tween intejaarticle interactions and external periodic 
potentials Jim 

Double-layer quantum Hall (2LQH) systems can show 
both types of behavior: they suraaort unusual fractionally 
charged topological excitations,BE3 and they apparently 



exhibit a CI transition to a state with broken transla- 
tional symmetry in theupresence of a sufhcently strong 
in-plane magnetic fieldJj It is interesting that in 2LQH 
systems, broken translational symmetry pan coexist with 
the hidden off-diagonal long-range ordenlil characteristic 
of quantum Hall states. The evidence for a CI transition 
in 2LQH systems is, however, indirect. Activation en- 
ergy measurements show that a sufficiently large in-plane 
magnetic field drives a transition between two different 
types of many-body ground states,t3 and the strength of 
the in-plane magnetic field and,±he size of the gaps are 
consistent with the CI scenario Jala But in order to estab- 
lish the nature of the competing quantum Hall ground 
states and the transition between them, additional types 
of experimental measurements are needed. Various types 
of experimental signatures of the CI scenario have been 
proposedp— including predictions of the form of the ep. 
ergy gapjij a predicted Kosterlitz-Thouless transition,E£l 
anisotropic transport in narrow samples,! 14 ! aiid the field 
dependence of the in-plane magnetization Here we 
propose that the CI transition in 2LQH systems could 
be studied directly by capacitance measurements, and 
by anisotropic transport produced by the dipole-density 
wave described below. 

This paper examines a novel effect of unequal layer 
densities on the CI transition in 2LQH systems. (In- 
terestingly, unequal layer densities can actually enhance 
the stability of the 2LQH state at a total filling factor 
i/>r=l.t3c3) It is found that when the layer densities are 
not equal, they go from being uniform in the commen- 
surate phase to becoming "rippled" in the incommensu- 
rate phase; this allows the CI transitiop_.to be detected 
by sensitive capacitance measurements. E3 The layer im- 
balance can be produced in two ways: most commonly 
by an external bias, but perhaps also spontaneously in a 
tilted sample with an unusually small capacitive char.gipgj 
energy and a sufficiently large interlayer tunneling£3'EHI 
We will focus here on bias-driven imbalance, since it is 
easier to achieve experimentally. In addition, capacitive 
techniques provide a quantitative measure of the inter- 
layer exchange and pseudospin stiffness in 2LQH systems. 
This is illustrated in the following sections by a Hartree- 
Fock calculation of the layer densities and capacitances 
for a 2LQH state at a total filling factor i/t=1. Similar 
effects should, in principle, occur at other filling factors, 
although i^T=l is probably most promising for experi- 
mental observation. 
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II. RIPPLED STATE 

Interlayer Coulomb interactions at low filling factors 
can stabilize 2LQH states when the layer spacing is com- 
parable-ito the separation between electrons within the 
layers.Eil Even at a total "integer" filling factor i*r=l, ex- 
periments indicate that the quantum Hall ground states 
are stabilized by Coulomb interactions and do not require 
interlayer tunneling for their existence. EJEil Further evi- 
dence of the rich variety of 2LQH states at Ut=1 comes 
from measuring the effects of an in-plane magnetic field, 
which induces a transition between two types of quantum 
Hall ground states£3 These effects have been discussed 
in terms of an unusual broken-symmetry quantum Hall 
ground states that-eshibit spontaneous interlayer (phase) 
coherence (SILC)BH 

At sufficiently small layer separation, a 2LQH,|Sjj-stem 
is an unusual quantum itinerant ferromagnet.ffluj The 
SILC 2LQH quantum ferromagnet exhibits a rich variety 
of ground stateSy-ohase transitions, and charged and neu- 
tral excitations.tSEia Murphy et al. investigated the effect 
of an in-plane magnetic field B» on 2LQH systems, and 
found evidence of a phase transition between two-compet- 
ing QH ground states at a critical value Bu — B C S3 These 
two ground states have been explained theoreticallyu by 
showing that application of a sufficiently strong parallel 
magnetic field B» > B c produces a soliton-lattice (SL) 
ground state in the incommensurate phase of the 2LQH 
system. Recent measurements of the interlayer tunneling 
conductivity in bilayer quantum Hall samples have prp=, 
vided dramatic evidence for interlayer phase coherence.Eil 



The first two terms of Eq n £2-4) constitute the 
Pokrovsky-Talapov (PT) modelETEKj with coefficients 
that, in the Hartree-Fock approximation (HFA), depend 
on m z according to 



t ee to^I^fe-^ 2 / 4 = toVM^e"^ 2 / 4 (2.2) 
p s ee p (l - ml) = 4z^tA>Po, 

where to is the tunneling-matrix element when ^1=^2; it 
is equal to half the symmetric-antisymmetric gap As^s. 
In the presence of a parallel magnetic field, 



*o^toexp(-Q^74) 



(2.3) 



which is a single-body effect. El The interlayer pseudospin 
stiffness when the layers are balanced is 
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in the HFA. Thc-iialue of p s will be reduced due to quan- 
tum fluctuationsE3c3 and finite-thickness effects. By ad- 
justing the front and back gate voltages of the sample, v\ 
and Vi may be varied (with i/p e^-jji + ^2 = 1), thereby 
allowing t and p s to be adjusted.EZI 

The th ird term (quadratic in m z ) in the energy den- 
sity [Eq. (f2j])1 is a cagacitive charging energy that favors 
equal layer densities n The capacitive energy U is given 
in terms of the electrostatic Hartree energy {£>i) and the 
intralayer (Eq) and interlayer (E\) exchange energies by 



U = D 1 -E +E 1 



(2.5) 



A. Effective Hamiltonian 

Formally, it is simplest to obtain the ground-state char- 
acteristices of the SILC 2LQH state from the energy per 
unit area within a gradient approximationEi in which the 
pseudospin m(r) is assumed to vary slowly on the scale 
of the magnetic length i. In doing so, it is convenient to 
specify the order parameter m(r) in terms of two quanti- 
ties: m z (r), the local difference in layer occupancies, and 
0(r), the projected angle of m in the xy plane measured 
with respect to the x axis. For constant m z —v\-vi and 
in-plane magnetic field B11, the energv_per unit area of 
the SILC 2LQH state has the fonrP*™ 



D x = 



£ = 
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\ v a m z], 



(2.1) 



where £ has been expressed in terms of 8 = 9 + Q • r, 
Q = z x B||27rd/0o, d is the interlayer spacing, and cj>o = 
h/e is the magnetic- flux quantum. Mean-field equations 
for 9 and m z are obtained by minimizing Eq. (2.1) with 
respect to those same quantities. 
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where the exchange integrals Ij have been evaluated in 
the HFA. Most treatments of the SILC 2LQH state have 
been for equal layer densities (m 2 =0), since there is a 
significant cost in capacitive charging energy to unbal- 
ance the layersJj However, an application of back and 
front gate voltages allows charge to be transferred from 
one layer to another, giving rise to a tusable nonzero 
value for m z in the 2LQH ground stateB3 The effects 
of charge-transfer imbalance were studied both theoret- 
ically and experimentally, and it was found that charge 
imbalanGe,ean actually increase the stability of the 2LQH 
state-IHIEaEd 

We shall estimate the numerical values of our re- 
sults for a hypothetical— ."typical" GaAs (m* rj 0.07m e , 
e r w 13) 2LQH samplejl3 with a total density nx=l-0 x 
10 11 cm -2 , a layer (midwell to midwell) separation 
d=20 nm, and a tunneling energy io=0.5 meV (AgAS = 
11.6 K). Such a sample would have I sw 12.6 nm, djl s» 
1.6, Tiuj c sa 6.9 meV for v T =l, and_e 2 /47re£ w 8.8 meV. 
In the HFA, p sa 0.03 meV, D x = 14 meV, and 
U = 7.4 meV. 
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The effect of the gate voltages is described by the last 
(V g ) term in the energy density in terms of effective filling 
factors Pf and v-q for the front (F) and back (B) gates: 



Q s /Q c = {it/2f/[r 1 K{i 1 % 



(2.6) 



The effective filling factors v a (where cv=F,B) are defined 
by the the electric fields E a produced by the front and 
back gates through Gauss' law E a = ei> a /2ir£ 2 e, where e 
(approximately 13eo for GaAs) is the dielectric constant 
appropriate to the 2LQH sample. 



B. Parallel magnetic field 

When Q 0, the pseudospin stiffness p s competes with 
the rotating Zeeman pseudofield t cos 9 to determine the 
spatial orie ntations of the pseudospins. Minimizing £ 
in Eq. (2.1) with respect to 9 gives the two-dimensional 
sine-Gordon equation 



£ 2 V 2 (9 



(2.7) 



where the width of the soliton is proportional to the 
length, 



£ = y/2n Pa /t = £o (1 



2\l/4 



(2.8) 



which for our hypothetical sample gives £o ~ 17 nm. 
The soliton width £ sets the scale for spatial variations of 
the pseudospin m(r); thus the condition for the validity 
of the gradient approximation is that £ be significantly 
larger than the magnetic length £. 

For sufficiently small Q, Eq. (2.1) is minimized by 
9(r) = 0. This is called the commensurate (C) phase, 
and in this phase the pseudospins align themselves with 
the rotating Zeeman pseudofield, so that 9 = — Q • r. 
However, above a critical value of B\\ corresponding to 

Q c = 4/(tt£) = (4/ir£) y/t/2irp s , it becomes energetically 
favorable to produce dislocation lines (solitons) . The soli- 
ton widths are of order £. Solitons proliferate rapidly 
for Q > Q c (incommensurate phase) because they re- 
pel each other only very (exponentially) weakly. The 
resulting array of solitons breaks the translational sym- 
metry of the 2LQH ground state by forming a SL. For 
large Q S> Q c , the rapidly varying tunneling phase factor 
causes the pseudospins to behave (nearly) as if t=0. In 
the HFA, the critical value of the in-plane field varies with 
the layer filling factors like Q c cx (1 — to 2 ,)^ 1 / 4 ; thus tun- 
ing the layer filling factors (m z ) via gate voltages allows 
the location (Q c ) of the CI transtion to be fine tuned.EJ 
The density of soliton lines in the incommensurate 
phase is proportional to the soliton wave vector Q s = 
2ir/L s , where L s is the spacing between solitons in the 
SL. The soliton density (proportional to Q s ) is calcu- 
lated as a function of the in-plane field (proportional to 
Q) via two equations involving an inJjswnediate parame- 
ter 77. The parameter r\ is defined byuO 



(2.9) 

and approaches 1 near the CI transition (Q s — > 0) and 
goes to deep in the incommensurate phase (Q s — > 00). 
Here -ST (77) is the complete elliptic integral of the first 
kindiil 
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FIG. 1. Soliton density (times 2-7r) vs in-plane magnetic 
field for balanced layers (m z —0, solid curve) and for the rip- 
pled state (m 2 o=0.5, dashed curve). The soliton density rises 
abruptly in the incommensurate phase, Q > Q c . There are 
no solitons in the commensurate phase, Q < Q c . 



When the layers are balanced (v\ =jl/; 
total energy with respect to Q s giveso' 

Q/Q c = E(t])/ij, 



minimizing the 



(2.10) 



whei fJH(ri) is the comp lete elliptic integral of the second 
kind.L 1 ] Equations ( |2.10| ) and ( 2.£ ) together determine the 
soliton density (or Q s ) as a function of the in-plane mag- 
netic field (or Q). Q s /Q c is plotted as a function of Q/Q c 
in Fig. 1, for the balanced casal3 (m z =0, solid curve) and 
for the rippled state (m z =0.5, dashed curve). Note the 
abrupt rise in soliton density for Q > Q c . The compres- 
sional stiffness K\ of the Sl£3 is proportional to the slope 
of the Q s versus Q curves in Fig. [TJ (see Ref. |l6| ): 



Ps 



dQ s 
dQ 



(2.11) 



As expected, Kx vanishes as Q — ► Q c 



C. Rippled layer imbalance 

Minimizing Eq. (|]l]) with respect to variations in 
gives 
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V a 



U+ [2tcos6- 4ir£ 2 p s (\'9 



Q) 2 ]/(l-m 2 ) 

(2.12) 



This equation determines (in the Hartree-Fock gradient 
approximation) the filling factor of each layer for the 
!/j = 1 SILC state. For sufficiently small in-plane mag- 
netic fields (Q < Q c ), 9=0 (the commensurate state), and 
the assumption that.m z is constant is self-consistent, pro- 
vided that t -C UsEi However, when Q > Q c , the quan- 
tum Hall ground state breaks translational invariance: 9 
is not spatially uniform and the soliton-lattice state is 
obtained^ When Q > Q c and v\ ^ i>2, Eq. ( |2.12 ) shows 
that m z also breaks translation invariance, and one ob- 
tains "rippled" layer densities. Thus uniform m z is not 
consistent with the broken translation symmetry of 9=0 
in the incommensurate phase, when T>\ ^ v^. The result- 
ing behavior of m z (r) is illustrated in Fig. [| 
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FIG. 2. Rippled layer imbalance m z (x) vs position for a 
soliton lattice with spacing L a /£=20 between soliton lines. 
We have used the parameters for the "typical" sample de- 
scribed in the text, and taken m z o=0.5. For these parameters, 
m z i=0.048. 



It is important to note that Eq. ( 2.12| ) is valid even 
in the incommensurate state of spatially varying 9, but 
only up to first order in t/U. Since t/U is small in the 
bilayer samples studied so far, this assumption is not too 
restrictive. From now on, we shall work only to lowest 
non-trivial order in t/U, and expand m z as described in 
the Appendix. Thus we shall not include the effect of the 
rippling of m z on the soliton width £, since this would 
produce only small corrections (i.e., of higher order in 
t/U) to t he results presented here. The soliton width £ in 
Eq. (2.8) is therefore computed using m z » m z o, where 
m Z Q is the layer imbalance in the absence of interlayer 
tunneling, as defined in the Appendix. The value of m z o 
depends on the gate voltage V g , but not on the "rippling" 



effect (which is of higher order in t/U). 

Because m z =v\-vi is associated with differences in 
layer electron densities, the rippling has the effect of asso- 
ciating an electric-dipole density with each soliton. When 
the solitons are separated (L s > £), the dipole- moment 
per unit length is (see the Appendix) 



Sp 



ed£o 2t 



m z0 



2itl 2 U 



(1 



L z0 



,1/4 



(4Q/Q c -3), (2.13) 



which has a value of about — O.lOe when to z o=0.5 and 
Q=Q C - Because the solitons have associated dipole mo- 
ments, the dominant interactions between solitons when 
the solitons are separated will be their dipole-dipole re- 
pulsion. In the limit that the solitons are well separated 
(L a 3> £,d), the interaction per unit length between two 
solitons separated by a distance x is 



V _ (5p/5y) 2 
L 2-ntx 2 '■ 



(2.14) 



so that the solitons repel each other with a force that 
falls off with distance as an inverse power (for v\ ^ U2), 
rather than exponentially (for v\ = z/2). Summing all the 
the dipole-dipole soliton interactions gives, in the ther- 
modynamic limit, 



V 



L T L n , 



tt {Sp/Sy) 2 



12 



1 (5p/8yf 
96tt 2 e 



(2.15) 



which is proportional to (t/U) 2 and is therefore small in 
magnitude. However, near the CI transition, the soli- 
tons are well separated and dipole interactions dominate 
the repulsions between the solitons, which has a strong 
effect on the compressional stiffness K\ of the SL. The 
relation between the wave vector Q and the parameter r\ 
is obtained by minimizing the total energy per unit area 
with respect to Q s at fixed Q (Ref. |l6|); when the lay- 
ers ar e balanced, Eq. ( 2.1 0| ) results, and Eqs. (2.9) and 
( |2.10D may be combined to obtain Q s as a fu nction of Q. 
When the layers are imbalanced, Eq. (2.1C) acquires an 
additional term due to the dipole interactions between 
solitons, 



Q/Qc = E( v )/ v - 



1 



d V 



PsQc dQ s L x Ly 



(2.16) 



C(Q S /Q c f 



where C ~ 0.14 for the hypothetical "typical" sam- 
ple with m z p=0.5 and Q=Q C - Near the CI transition 
Eq. (|2~1~6|) gives 



Qs/Qc w ViQ/Qc-i)/C, 
K x / Pa a 2y/C(Q/Q c -l). 



(2.17) 



The corresponding formula for the balanced (^1=^2) case 
are quite different: Q s ~ — l/ln(Q — Q c ) and K\ ~ 
(Q ~ Qc), up to logarithmic corrections. Nonetheless, 
K\ (and, by definition, Q s ) vanishes at the CI transition 
which has important consequences for the capacitance 
near the CI transition. 



4 



D. Topological charge 

We also note that when the layers are imbalanced, 
the soliton lines acquire a charge density when they are 
tilted or sheared. In the lowest Landau level, the pseu- 
dospin textures with topological charge possess an elec- 
tric charge. The electric-charge density Sp(r) is just the 
topological-charge density (i.e., the Pontryagian index 
densityM) times —evr, 



Sp = — — e^m • (dim x 9,-m), 

07T 



(2.18) 



where is the totally antisymmetric tensor of rank 2, 
and we sum over i and j over the values 1 and 2 [x 
and y). The tot al fil ling factor is defined as Vt=V\-\-v 2 . 
Expression Eq. (2.18) in terms of m z and 9 using 



rax cos ( 



h = VI- 



(2.19) 



Pxx ~ Pyy: Q ^ Qc-t 

Pxx ^ Pyy-; Q Qc-t 



(2.24) 



Thus a disparity between the longitudinal resistivities 
parallel and perpendicular to the in-plane magnetic field 
indicates the presence of an incommensurate soliton- 
lattice phase in the imbalanced system. 

We note that anisotropic transport in the quantum 
Hall regime has been found in a very different context, 
in very high-mobility, skiff le-layer samples at high half- 
integer filling factors.t3t3 In that case, the anisotropic 
transport is taken to indicate the presence of a spon- 
taneously striped anisotropic charge density of the quan- 
tum Hall ground stateO Here, the stripes are not sponta- 
neous (for t <C U) , but are produced and oriented by the 
in-plane magnetic field. Nonetheless, the stripes (dipolar 
soliton lines) found here should also produce anisotropic 
transport. 



gives 



Sp(r) = ^(d x 6d y m z 



d x m z d y 6 



(2.20) 



4tt 



(d x 8 - Q)d y m z - d x m z d y 



where we have taken Q to be along the x direction. If 
we rotate the soliton lines by taking 



h(x) -> 9 (ax + (3y), 



(2.21) 



and also transform the rippled layer imbalance according 
to 



5m z i(x) — ► 6m z i(ax + f3y), 
then the associated charge density is 

CUT 

Sp = — — — Qd y Sm z i(ax + (3y), 

47T 



(2.22) 



(2.23) 



which associates a spatially varying charge density pro- 
portional to (3t/U to the tilted soliton lines. The inte- 
grated charge_of the soliton lines remains zero, so the 
shear stiffnesses may not be greatly affected. 



E. Anisotropic transport 

Application of a sufficiently strong in-plane magnetic 
field produces soliton lines that are parallel to the in- 
plane magnetic field. When the layers are imbalanced, 
these soliton lines possess dipole-charge densities. The 
resulting electric fields associated with the dipole-charge 
densities will make the conductivity anisotropic in the 
incommensurate phase. Transport parallel to the soliton 
lines (along the direction of the in-plane magnetic field, 
sa Y — y) is expected to be easier than in the direction 
perpendicular (along x) to the soliton lines: i.e., when 
the layers are imbalanced, we expect 



III. INTERLAYER CAPACITANCE 

Capacitance measurements offer the possibility of 
directly probing the ground-state properties of two- 
dimensional electron systems, such as the thermody- 
namic compressibility. Although the differential gate 
capacitance, which measures how the charge on a gate 
changes with respect to changes in the gate voltage, is 
slightly affected by the compressibility of the electron gas 
in the occupied layer that is nearest the gate, it is almost 
entirely dominated by the large gate-to-layer distance of 
the device. A far more sensitive measurement of the elec- 
tronic compressibility is provided by the. Eisenstein ratio 
Re which is an interlayer capacitance,o~E3 



Re = SE 12 /SE t 



gate j 



(3.1) 



where E\ 2 is the electric field that exists between the 
layers, and -E ga tc is the electric field between the gate 
and the nearest layer. Classically, conduction electrons 
in the layers should completely screen the electric fields 
produced by the gates, so that Ei 2 =Q and Re=0- Indeed, 
this result is approached when the layers are sufficently 
far apart (beyond several hundred A). But, due to their 
finite density of states, the effectively two-dimensional 
electron layers cannot completely screen the gate electric 
fields, so E 12 is nonzero when the gates are unbalanced. 

The Eisenstein ratio has been measured experimen- 
tally in two-dimensional electron gas (2DEG) systems, 
and used to determine the negative compressibility of the 
low-density 2DEG^3 it haSjaJsa been calculated theoraU 
ically at zero magnetic fieldEjO and in 2LQH systems.E£l 
The calculation of the Eisenstein ratio is discussed in 
some detail in Refs. |35|^7|; here we shall briefly outline 
only the key steps. 

It is covenient to separate the total energy per unit area 
of the 2LQH system into an electrostatic part (i.e., the 
integrated electrostatic energy density eEf 2 d/2 between 
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the layers) and a many-body part (e), which would be 
the energy per area for a sjptem with neutralizing charge 
backgrounds in each layer The chemical potential rel- 
ative to the bottom of quantum well j is given byEJ 



Hi = d{e)/drii 



(3.2) 



where rij — Vijlul 2 is the areal number density in layer i. 
Re can be expressed in terms of the interlayer separation 
d and the effective electronic lengths Sy, defined as 



_ ^(Hm_ _ I d 2 / 27r^ 2 (e) \ 
~ e 2 drij ~ 2 di/jdut \e 2 /4ire£ J 



(3.3) 



from which it follows that Sji = sy. In the absence of 
interlayer interactions (the case considered in Ref. 35), 
S12 — s 2i — and the length su is inversely proportional 
to the electronic compressibility Ki in layer i (Ref. |36| ): 

(3.4) 



e 2 n 2 Ki ' 



If the front-gate and back-gate voltages are varied si- 
multaneously so that the total layer density is kept fixed 
(i/i+V2=l), then 



Re — SE 12 /SE gate = 



si + s 2 



d + si 



S2 



where 



Sl = Su - Sl2, 



S2 = S22 - S21. 



(3.5) 



(3.6) 



If only one of the gate voltages is kept fixed 
numerator of Eq 



then the 

5j) is equal to either si or sa, in- 



stead of their sum.KZTAt high densities, Sj (and thus Re) 
are positive, but at suffiently low densities, they become 
negative, resulting in the negative values of Re that have 
been measured experimentally.^ 

The Eisenstein ratio has been calculated for a vt = 1 
2LQH state,E3 although without a parallel magnetic field. 
It was found that, although the lengths s, were negative, 
the criteiipn for stability against abrupt interlayer charge 
transfer E3 



d + si + s 2 > 0, 



(3.7) 



is still satisfied. It follows from Eqs. ( |3.5| ) and (3/7) that 
when Re diverges, it signals an interlayer charge-transfer 
instability. 

In the next two subsections, we calculate the contribu- 
tion Rei of the PT Hamiltonian to the Eisenstein ratio. 
Rei contains the dependence of Re on the tunneling t, 
pseudospin stiffness p s , parallel magnetic field Q, and 
layer imbalance m Z Q. It is convenient to separate Rei 
into two parts: one that is formally divergent at Q=Q C 
when the layers are unbalanced (m z 0), and one that 
does not diverge, but which nevertheless exhibits a non- 
trivial dependence on Q and m Z Q. 



A. Divergent contribution 

We shall first consider only the contributions to the 
electronic lengths Si that become negative and divergent 
at the CI transition. We therefore focus on how changes 
in the layer imbalance m z allow one to cross through the 
CI transition, due to the m z dependence of the critical 



parallel ma gnet ic field (or of Q c ), as discussed in Sec. II B 
and in Ref. |l6|. We therefore consid er o nly the PT part 



of the the energy per area in Eq. (2.1), and then only 
that part of the PT energy that depends on m z through 
Q /Qc- (There are other terms that contribute to s,, but 
they do not give rise to a divergence at the CI transition, 
at least for t/U <C 1, which is the usual situation in 
most samples.) Because we are focusing on the effects 
of passing throught the CI transition, rather than the 
effects of changing the total filling factor, we restrict our 
attention to the case of fixed total filling factor vt- With 
these restrictions, the effect of taking a derivative of a 
function of Q/Q c with respect to the layer filling factor 
Vi produces the equivalence 



d dm z dQc Q d 

dvi dvi dm z Q c dQ 



(-1) 



d 



2 l-mfdQ- (3 - 8) 



Combining Eqs. (3.3) and (3.8) with the magnetization 
calculations of Ref. [16| gives 



si + s 2 



1 — mi 



{Qtf 



dQ s 



e 2 /4TTe£ dQ 



(3.9) 



where we have kept only the contribution that diverges at 
the CI transition. The divergence occurs because nonzero 
m z allows small changes in the layer filling factors to 
tune the system through the CI transition. The divergent 
part of the electronic lengths Sj are thus proportional 
to the in-plane differential magnetic susceptibility that 
was calculated in Ref. |l^, and which diverges at the CI 
transition. 

Unfortunately, the divergence in the electronic lengths 
Si at the CI transition is weak, proportional to 
\j \/Q — Qc, with a small prefactor. For the hypothet- 
i cal s ample con sidered in this paper and m z o=0.5, Eqs. 
~Ti| and (lO) give 



(I 



si + s 2 

I 



0.005 



VQ/Qc 



l 



(3.10) 



5 to give a divergent 
d + si + S2 = 0), at 



which requires (Q/Q c — 1) ~ 10 
value of the Eisenstein ratio (i.e. 
least for a tunneling-matrix element < =0.5 meV. The 
required nearness to the CI transition (i.e., the smallncss 
of Q/Qc — 1) is inversely proportional to to, so making to 
smaller might be of some help. However, it may be that 
in practice, sample disorder smears out the CI transition 
well before the divergence in Re can be approached. 



G 



B. Nondivergent contributions 



We saw in Sec. Ill A that, although there is a contri- 
bution to the interlayer capacitance (the Eisenstein ratio 
Re) which is formally divergent at the CI transition, it 
may be difficult in practice to tune close enough to the 
transition to observe the divergence. It may be that the 
nondivergent contribution of the PT contribution to Re 
is more easily detected. We calculate this contribution 
below for small 2t /U. 

For simplicity, we focus on the case of fixed total fill- 
ing f acto r. Then it follows from the definition of Re in 
Eq. (3.1), from Gauss' l aw, fr om the definition of V g in 



Eq. (2.6), and from Eq. ( 2.12 ), that we may express Re 
as 



Ri 



1 



Di d(r 



U dm z0 
R E0 = 1- Dx/U, 

-Di d(m zl ) 



Reo + Rei 



(3.11) 



Rei 



U dm. 



2to 
U 



(cos 6» ) 



(l-m*o)-3/= 



-^((V0 o -Q) 2 



The largest contribution to the Eisenstein ratio is Reo, 
and it was this quantity that was calculated in Ref. |3(]. 
For the hypothetical sample parameters in the text, 
Reo = —0.9, independent of m z o and Q. Re = Reo 
when the interlayer-tunneling amplitude to is very small, 
or when the parallel field is large (Q substantially larger 
than Q c ). 

Rei contains all the dependence of Re on the in-plane 
magnetic field and layer imbalance. If t/U is made as 
large as possible, then by measuring the Q and m z o de- 
pendence of Re, it may be possible to measure Rei- In 
the commensurate phase (Q < Q c ), 9o = 0, so 



Rei 



Di 2t /U 
U y/l 



'zO 



1 /4_Q 

(!-"&)) \kQc 



(3.12) 



In the absence of an in-plane magnetic field (Q—0), 
D 1 2t /U 



R B1 (Q = 0) = - 



U 



1 



3/2 
l z0 



-0.39, (3.13) 



for the hypothetical sample parameters described in the 
text, with m,zo=0.5, so that in this case Rei is almost 
half the size of Reo- As Q increases, Rei decreases in 
magnitude and, provided that m z o < \Jl — (7r/4) 2 w 
0.62, passes through zero and becomes positive near the 
CI transition. 

At the CI transition, this nondivergent part of Rei is 

2 



R E l{Q = Qc) = 



Di 2t /U 



m 



zO 



1 



(1 



4) 



(3.14) 



which for the sample parameters described in the text has 
a value of Rei(Q = Q c )=0.16 for m z o=0 and Rei(Q = 
Q c )=0.08 for m z o=0.5. In the incommensurate phase, 
Rei drops rapidly to zero due to rapid spatial variations 
in #o (x)Jl3 By measuring how the in-plane magnetic field 
and layer imbalance affect the interlayer capacitance Re, 
it may be possible to estimate the pseudospin stiffness 
p s in the commensurate phase, and also to detect the 
incommensurate phase, as signaled by a rapid decrease 
in the sensitivity of Re to Q and m zQ . 



IV. CONCLUSIONS 

It has been shown that the incommensurate phase of a 
bilayer quantum Hall state has a "rippled" dipole charge 
density whenever the layers are unbalanced. The rip- 
pling arises because the layer imbalance m z and the in- 
terlayer phase 9 are coupled through the m z dependence 
of the effective tunneling energy t and pseudospin stiff- 
ness p s in t he Pokrovsky-Talapov part of the total energy 
in Eq. ( |2.1[ ). This coupling between the layer imbalance 
and interlayer phase produces the rippled state when the 
translational symmetry of the phase is broken in the in- 
commensurate phase. The rippled layer imbalance was 
calculated within the Hartree-Fock gradient approxima- 
tion, and is illustrated in Fig. ^. The details of the calcu- 
lation are given in the Appendix. We focused on the limit 
where the interlayer-tunneling energy is smaller than the 
charging energy of the bilayer (t <C U), which is the case 
for all bilayer samples which have been studied experi- 
mentally so far. 

Because solitons have an associated dipole-moment per 
unit length (which we estimated in the Appendix) in the 
rippled state, well-separated soliton lines experience a 
power-law (inverse cube) repulsive-force per unit length, 
rather than the much weaker exponentially-decaying re- 
pulsion between solitons found in the balanced case. This 
has a strong effect on how the density of solitons depends 
on the in-plane magnetic field near the CI transition, as 
illustrated in Fig. [l| 

The fact that solitons have an electric dipole moment 
in the rippled state makes it likely that they produce 
anisotropic transport. Transport parallel to the soliton 
lines (±y) is likely to be easier than transport perpen- 
dicular to the soliton lines (±x). We expect the ratio 
Pxx/Pyy to increase when the rippled state is entered. 
This could provide an experimental signature for the in- 
commensurate state. 

We calculated the interlayer capacitance, specifically 
the Eisenstein ratio Re, which is a sensitive measure of 
the electronic compressibility and of interlayer electronic 
correlations. When the layers are unbalanced, there is 
a contribution to Re which diverges at the CI transi- 
tion. However, observing this contribution is likely to be 
problematic because it requires that the CI transition be 
very sharp (unsmeared by disorder). We also calculated 
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a contribution to Re which does not diverge at the CI 
transition, but which nevertheless offers a good possibil- 
ity of experimental detection. We calculated the in-plane 
magnetic- field dependence of Re, along with its depen- 
dence on layer imbalance, especially in the commensurate 
(Q < Q c ) phase, up to the CI transition [Q—Q c ) } and dis- 
cussed its behavior (rapid decline) in the incommensurate 
(Q > Q c ) phase. By measuring Re, the pseudospin stiff- 
ness could be estimated and the incommensurate phase 
detected. 

The existence of "rippled" layer densities in the in- 
commensurate phase of unbalanced 2LQH systems is 
expected to be valid beyond the HFGA, although the 
size of the density, variations will be reduced by quan- 
tum fluctuationsBU'Eij As long as the layer densities are 
"rippled" , they will make an anomalous contribution to 
the capacitance.. Alt hough fluctuation effects change the 
sizes of p s and t,E3E3 the basic physics of producing "rip- 
pled" layer densities still holds beyond the HFA. In par- 
ticular, p s and t will still change with the layer imbal- 
ance m z , and this dependence on m z will produce a cou- 
pling between m z and the pseudospin angle 8, leading 
to nonuniform m z in the incommensurate phase when 

Several remarks about the observability of the capaci- 
tive effects found here are in order, especially the effects 
of including finite temperature and disorder. In practice, 
both finite temperature and disorder limit the minimum 
effective value of (Q/Q c — 1) that can be obtained. These 
are important topics which deserve further study; only 
some preliminary considerations are discussed here. 

Ignoring the effects of disorder for the moment, it is im- 
portant to note that the soliton lattice exists only at suf- 
ficiently low temperatures. The soliton lattice supports 
a finite-temperature Kosterlitz-Thouless (KT) transition 
due to dislocation-mediated melting of the lattice of soli- 
ton linesE'ES The KT temperature for melting the soliton 
lattice is roughly /cbTa't ~ (tt /2)\/K\Ki, where K\ is 
the longitudinal stiffness ,-apd K 2 is the transverse stiff- 
ness of the soliton lattice £j Because K\ — > as Q — > Q c , 
the KT temperature drops as the CI transition is ap- 
proached. This sets a limit to how close one can get to 
the CI transition, which may be estimated in the case 
of finite l ayer i mbalance by usiug Ki — 2y/ C(Q/Q C — 1) 
[see Eq. (2.17)] and K 2 = p s -t3 From these considera- 
tions, the requirement that T < Tkt gives 



(Q/Qc 



l) > 



1 (2k B T 



7T p s 



2.4 



T 

Tk 



(4.1) 



For 



for the sample parameters used in the text. 
T=100 mK, this yields (Q/Q c - 1) > 2.4 x 1CT 4 . 

The smearing of the CI transition due to disorder is 
more problematic. Even in capacitance experiments de- 
signed to measure the (in) compressibility of the fractional 
quantum Hall (FQH) state (which would, in the absence 
of disorder, lead to a divergent Kg at odd-denominator 
filling factors), only finite changes in Re are found at the 



FQH filling factors, due to the effects of disorder.il It is 
to be expected that disorder may eliminate any abrupt 
features in the interlayer capacitance Re in this case also. 
Strictly speaking, the long-range order ofr-the SL is de- 
stroyed by any finite amount of disorder ;E3 presumably 
the SL has only a finite correlation length < oo due 
to disorder. Roughly speaking, the maximum spacing 
between solitons (L s ) will be limited to L s < Cd> so that 



(Q/Qc 



1) > C 



ZEll 

2 (d 



(4.2) 



In practice, such disorder could arise from small varia- 
tions in the local tunneling ampitudes t or the spin stiff- 
ness p s due to minute variations in the interlayer barrier 
thickness and/or the layer separation. More theoretical 
work needs to be done to determine the limits imposed 
by disorder, but the issue of the observability of the CI 
transition in capacitance measurements must be settled 
experimentally. 
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APPENDIX: DIPOLE MOMENT 

In this appendix, we estimate the dipole-moment per 
unit length associated with solitons in the incommensu- 
rate pha se of an imbalanced 2LQH system. Expanding 
Eq. (2.1) in powers of (t/U) gives 



m z (r) ps m z0 + m zl (r), 

where 

m z0 = Vg/U 

is the f=0 result for the layer imblance, and 
2t m z0 



m z i(r) 



U(l-m 2 zQ ) 



r(v0 o - Q) 



cos 6*o 



(Al) 
(A2) 

, (A3) 
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to first order in t/U. Here 6o(r) is the soliton-line solution 
to the PT model for the val ues of t and p s corresponding 
to setting m z —m Z Q in Eq. (|2.2|). 

For a single soliton, the lowest-order solution for 9 ist£l 



6q(x) = 4 arctanexp(:r/£), 



(A4) 



where we have taken Q to lie along the x direction. This 
gives 



1 — cos6>o(a;) 



£d x 8 Q (x) 



cosh 

2 



(A5) 



cosh(a;/£) 
It is convenient to express m z i(r) as 
m zl (r) = m zl + Sm z i(x), 

where 

2t m z0 



U (1 - m 2 z0 ) 



(Q 2 e i) 



(A6) 



(A7) 



is the value of m z i in the commensurate (#o = 0) 
phase. For the hypothetical sample described in the text, 
fh z i ~ 0.048 for m z o=0.5 and Q « Q c . We associate the 
spatially dependent part of m z (see Fig.i ||) with the soli- 
ton line: 



Sm z i(x) 



2t m z0 
U(l-m 2 z0 ) 



4QC 



cosh(:c/£) cosh 2 (x/£) 



(A8) 



to 




The areal number density of layer j is given by rij = 
fj/(2ir£ 2 ). Therefore, the dipole-moment per unit area 
is 



P 



L X Ly 



ed 
2^ 



rn. 



(A9) 



and we associate an dipole-moment per unit length 
ed 



Sp 
5y 



2tt£ 2 



5m z i [x)dx 

-oo 

ed£ 8t m z0 



(A10) 



m z0 



2< U (l-m 2 )i/4 



(4Q/Q c -3). 



For the "typical" sample described in Sec. II A at layer 
imbalance m Z Q = 1/2 and Q=Q C , 



Sp/5y w -O.lOe, 



(All) 



where — e is the electric charge of an electron. 

When the soliton lines do not overlap (Q sufficiently 
near Q c ), then 6q(x) is very nearly a periodic superposi- 
tion of single-soliton solutions, spaced apart by L s : 

§ (r) w 4^arctanexp[(a; - jL a )/g\. (A12) 



The dipolar interaction-energy per unit length between 
two parallel soliton lines separated by a distance x is 



V 2 (x) (Sp/Syf 



L y 2Trex 2 



(A13) 



Thus the total dipole-interaction energy per unit area is 
V N a 



(A14) 



3=1 

7T (Sp/Sy) 2 
12 eU 



1 (Sp/Sy) 2 



96tt 2 



10.0 



where we have used the fact that the number of solitons 
is iV s = L z /L s . 

The relation between the wave vector Q and the pa- 
rameter 77 is obtained by minimizing theJjptal energy per 
unit area with respect to Q s at fixed Q;E3 when the lay- 
ers are balanced, Eq. ( 2.10| ) results, and Eqs. (2.9) and 
( |2.10 ) may be combined to obtain Q s as a fu nction of Q. 
When the layers are imbalanced, Eq. (2.1C) acquires an 
additional term due to the dipole interactions between 
solitons, 



FIG. 3. Spatial dependence of the layer imbalance 
8m z i(x) associated with a single soliton, vs position for the 
hypothetical "typical" sample parameters described in the 
text. 



Q/Q c = E( v )/ V 



d V 



PsQc dQ s L x Ly 
E^/rj + CiQs/Qc) 2 , 



(A15) 



where 







8tt V 



1 fSp/Sy\ e 2 /Airel 



QJ, 



(A16) 



and C ~ 0.14 for the hypothetical "typical" sample with 
"izO—0. 5 and Q=Q C . Because the interaction between 
separated solitons is an inverse-power law (when unbal- 
anced) rather than exponetially decaying function (when 
balanced), Q s and K\ are proportional to y/Q — Q c near 
the CI transition. 
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